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STAR FORMATION AND TURBULENT DISSIPATION IN MODELS OF 

DISK GALAXY EVOLUTION 

V. Avila- Reese 1 , C. Firmani ' 2 , and E. Vazquez-Semadeni 1 
RESUMEN 

La tasa de disipacion de energfa cinetica en el MIE turbulento dc las galaxias de disco es un ingrdicntc clave 
para modelos de evolution galactica, ya que determina la efiticncia de la retroalimentacion de la formation 
estelar (FE) de gran escala. Con ayuda de simulaciones MHD encontramos que el MIE disipa eficientemente la 
cnergi'a turbulenta inyectada por fuentcs dc naturalcza cstclar. Asi, el proceso de FE puede ser auto-regulado 
solo a nivel del MIE del disco. La aplicacion del mccanismo de FE auto-regulada en modelos de evolution 
galactica, donde los discos se forman dentro de halos de CDM en crecimiento, permite predecir la historia de 
FE de las galaxias de disco, incluyendo a la Via Lactea y la vecindad solar, asf como la contribution de toda 
la poblacion de discos a la historia de FE cosmica. Los resultados son alentadores. 

ABSTRACT 

The kinetic energy dissipation rate in the turbulent ISM of disk galaxies is a key ingredient in galaxy evolution 
models since it determines the effectiveness of large-scale star formation (SF) feedback. Using MHD simulations, 
we find that the ISM dissipates efficiently the turbulent kinetic energy injected by sources of stellar nature. 
Thus, the SF process may be self-regulated by an energy balance only at the level of the disk ISM. The use 
of the self-regulation SF mechanism in galaxy evolutionary models, where disks form inside growing CDM 
halos, allows to predict the SF history of disk galaxies, including the Milky Way and the solar neighborhood, 
as well as the contribution of the whole population of disk galaxies to the cosmic SF history. The results are 
encouraging. 
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I. INTRODUCTION 

Star formation (SF) and feedback are key pro- 
cesses in the evolution of galaxies, but probably they 
are the worst understood of all. The large-scale 
SF cycle in normal galaxies is believed to be self- 
regulated by a balance between the energy injection 
due to SF (mainly SNe) and dissipation. Two main 
approaches have been used to describe the SF self- 
regulation in models of galaxy formation and evo- 
lution: (a) the halo cooling flow-feedback approach 
(White & Frenk 1991), (b) the disk turbulent ISM 
approach (Firmani & Tutukov 1992; 1994; Dopita & 
Ryder 1994; Wang & Silk 1994). 

According to the former, the cool gas is reheated 
by the "galaxy" SF feedback and driven back to the 
intrahalo medium until it again cools radiatively and 
collapses into the disk. This approach has been used 
in the semi-analytical models of galaxy formation 
(e.g., Kauffmann, White & Guidcrdoni 1993; Cole 
ct al. 1994). The reheating rate depends on the halo 
circular velocity V c : M r h oc M s /V", where M* is 
the SF rate (SFR) and a > 2. Thus, the galaxy 
SFR, gas fraction and luminosity depend on V c . In 
these models, the disk ISM is virtually ignored and 
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the SN energy injection is assumed to be as efficient 
as to reheat the cold gas up to the virial tempera- 
ture of the halo. A drawback of the model is that it 
predicts hot X-ray halos around disk galaxies much 
more luminous than observed (Benson et al. 2000). 

According to approach (b), SF is triggered by 
gravitational instabilities (Toomre criterion) and 
self-regulated by a balance between energy injection 
and dissipation in the turbulent ISM in the direc- 
tion perpendicular to the disk plane: ^sn^sn^* + 
^E,accr = (£ fl 7^)/(2£ d ), where j S n and e S N are the 
kinetic energy injection efficiency of the SN into the 
gas and the SN energy generated per gram of gas 
transformed into stars, respectively, £* is the surface 
SFR, TiE,accr is the surface mass accretion rate, and 
E ff and v g are the gas surface density and velocity 
dispersion, respectively. The key parameter in this 
equation is the dissipation time td- Now, the ISM is 
a turbulent, non-isothermal, multi-temperature flow. 
It is known that in the denser regions of the ISM, 
where stars form, the cooling time is short with 
respect to the dynamical times (Spitzer & Saved- 
off 1950; Elmegreen 1993; Vazquez-Semadeni et al. 
1995, 1996). The main effect of this short cooling 
time is to cause the medium to behave essentially as 
a polytrope, with the pressure P scaling with den- 



1 



2 



AVILA-REESE, FIRMANI & VAZQUEZ-SEMADENI 



sity p as P cx p 1 , where 7 is an effective polytropic 
exponent that determines the compressibility of the 
medium (Vazquez-Semadeni et al. 1996; Scalo et al. 
1998; Spaans & Silk 2000; Jenkins & Tripp 2001). 
Given the medium's compressibility, the star forma- 
tion rate is then determined by the mean turbulent 
kinetic energy (which produces the compressions), 
which in turn is limited by the kinetic energy dis- 
sipation rate. Therefore, the key parameter is the 
turbulent dissipation time, and a crucial question is: 
how efficiently does the disk ISM dissipate the 
injected by stellar sources?. 

2. TURBULENT DISSIPATION IN ISM 
SIMULATIONS 

To explore the ISM dissipative properties, Avila- 
Reese & Vazquez-Semadeni (2001) used numerical 
MHD simulations of turbulent compressible fluids 
resembling the ISM of disk galaxies at the 1 kpc 
scale (Vazquez-Semadeni et al. 1995; Passot et al. 
1995). The MHD and the internal energy conserva- 
tion equations were solved in two dimensions 3 with 
model terms for radiative cooling and background 
heating. A pseudospectral method with hypervis- 
cocity and periodic boundary conditions was used. 
The way in which Ek is injected to the fluid is cru- 
cial: it should resemble "stellar" sources like SN rem- 
nants and expanding HII regions. That is, the en- 
ergy sources are discrete, generally small in compar- 
ison with the characteristic scales of the global ISM, 
and with filling factors typically <C 1. The forcing in 
the simulations was generated by "winds" : the gas 
around the SF regions received an acceleration af, 
directed radially away, thus producing an evolving 
velocity profile with characteristic forcing velocity Uf 
and radius If. The characteristic energy injected is: 
ekoc Ufaflf 2 A£ob, where, AioB is the source life- 
time (assumed constant and =6.8 Myr). The two 
independent input parameters that define the local 
energy injection process are If and ek (or Uf); at is 
already a response of the flow. There is a third ini- 
tial parameter related to the global energy injection: 
the surface SFR, E*. 

We performed a sizeable number of simulations 
to explore how the dissipative properties of the tur- 
bulent ISM depend on the forcing parameters: If, ek, 
and E*. In these simulations, the sources are placed 
randomly in space with a fixed probability, which 
determines the E*. 

The simulations show that the kinetic energy in- 
jection and dissipation rates are always very close 

3 We argue that dissipation in 2-D is representative of that 
in 3-D as long as it is dominated by shocks rather than by a 
turbulent cascade. 



to each other, suggesting that most of the dissipa- 
tion occurs at or near the sources, where shocks 
are common, and that, for practical purposes, the 
injection and dissipation times are equal. In gen- 
eral, a flow with physical parameters close to the 
Galactic disk at the solar neighborhood dissipates 
its turbulent energy rapidly, in approximately 15- 
20 Myr. A more realistic simulation, where the SF 
is density directed, gives td« 18 Myr (in this simu- 
lation, the input sources tend to cluster due to su- 
pershell formation). We also showed that, in terms 
of measurable properties of the ISM, td can be ex- 
pressed as: td> (Eg) u rms 2 /(ck E*), where (E g ) is 
the average gas surface density, and the inequality 
arises because some overlapping of the forcing re- 
gions may occur. Since E*oc E g " (Schmidt law), 
tdoc E g (" -1 )(r). Thus, td is expected to decrease 
with Galactocentric radius. 

Given the discreteness and low filling factor of the 
energy sources, the turbulence "propagates" away 
from them into the general ISM, reaching out to far- 
ther distances for lower td- The velocity dispersion, 
u rms , of the residual turbulence is found to decay 
roughly inversely with distance to the source. It was 
also found that the dominant parameter determin- 
ing td is If, while the energy per source and the 
SFR arc not very relevant (t d oclf - 84 /(c k E*) - 12 ). 
Thus, the production of very large active turbulent 
zones (able to inject energy out from the disk into 
the halo, for example), which occurs at large td, re- 
quires large values of the source characteristic size If, 
i.e., that the forcing regions are themselves large (su- 
pershells). These supershells could arise as a conse- 
quence of strong clustering and self-propagating SF. 
This might be the rule for starburst galaxies (non- 
stationary SF), but generally the exception in nor- 
mal disk galaxies. 

Our results, if applicable in the vertical direc- 
tion, imply that the turbulent Ek produced in the 
disk plane is mostly dissipated inside the disk ISM. 
The turbulent motions will propagate up to distances 
close to the observed HI disk semi-thickness, and will 
not be able to reheat and drive back the gas into the 
halo. The shock-dissipated turbulent energy trans- 
forms probably into thermal energy, which could be 
a source of high-energy photons able to ionize the 
low-density extraplanar medium (up to a few kpes 
from the plane). 

3. THE STAR FORMATION HISTORY 
A scheme of SF triggered by disk gravitational in- 
stabilities and self-regulated by the SN energy input 
and the turbulent energy dissipation was applied by 
Firmani, Hernandez & Gallagher (1996) to galaxy 
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Fig. 1. a). SF history in the solar neighborhood: pre- 
dicted by our model (thick line) and infered from ob- 
servations (jointed circles). Dashed line is for a MW 
model with a soft core in the halo. Dot-dashed line is 
the predicted gas infall rate at Rq . The vertical and 
horizontal shaded boxes are observational estimates for 
the present-day SFR and for the age of the solar neigh- 
borhood, respectively, b). Predicted evolution of the star 
(solid line) and gas (dashed line) surface densities at the 
solar neighborhood. 

evolution models. A Schmidt law with index < 2 
along a major portion of the disks was predicted. 
For closed models, even with the lowest disk surface 
densities, the SF timescales result shorter than 5 Gyr 
and the integral (B — V) colors redder than 0.9. Gas 
infall is necessary. This infall may be provided by the 
cosmological mass aggregation predicted in a hierar- 
chical CDM scenario, where disks form inside-out. 

Self-consistent models of disk galaxy formation 
and evolution within growing CDM halos, which use 
the disk self-regulated SF approach, indeed predict 
realistic integral and radial color indices, gas frac- 
tions, surface brightness profiles as well as realis- 
tic dynamical properties and correlations among the 
galaxy parameters (see Firmani & Avila-Rccse, this 
volume). In these models, the turbulent dissipa- 
tion time was calculated as td(r) = a(2f2(r))/fv(r) 2 , 
where fl(r) and n(r) are the angular velocity and 
cpicyclic frequency at radius r, and a is a constant 
close to unity (Firmani et al. 1996). The formula 
above was derived from a simple dynamical anal- 
ysis. The MHD simulations of the turbulent ISM 
presented in §2 offer the posibility to fix the value 
of a. Normalizing to the Galactic disk at the solar 
radius, we find that a ~ 1.5 — 2. The dependence of 



td on the radius in the formula above agrees with an 
inference from the simulations of such a dependence 
(see above). 

The SF histories of most of the model galaxies 
have a broad maximum at z ~ 1 — 2.5, and then 
the SFR decays by factors of 1.5-4 until z = 
(Avila-Reese & Firmani 2001). The SF history is 
driven by the halo angular momentum and by the 
halo mass aggregation history; in contrast with the 
halo cooling-flow model, there is not a significative 
dependence on mass or circular velocity. The an- 
gular momentum determines the disk surface den- 
sity. The higher the density, the more efficient is 
the SFR. The mass aggregation history determines 
the gas infall rate; for mass assembly histories of 
the halo more extended in time, the late gas accre- 
tion rate is higher, and therefore, the SFR is higher 
at z f« 0. Thus, HSB and/or red galaxies have a 
maximum SFR at earlier epochs and a faster decline 
towards the present than LSB and/or blue galax- 
ies. The latter may even have an increasing SFR at 
3 = 0. When comparing the SFR of high-redshift 
disk galaxy populations with the present one, it is 
important to take into account these differences of 
the SF history on galaxy surface brightness and color 
and the selection effects. 

3.1. Milky Way (MW) and the solar neighborhood 

The global and local SF histories of a MW model 
can be predicted. For a ACDM cosmology, the 
MW model was constructed using a virial mass of 

2.8 10 12 M Q , a baryon fraction of 0.02, the average 
halo mass aggregation history, and a spin parameter 
A = 0.02. The surface brightness profiles and scale 
radii (bands B and K), the global gas and stellar 
mass fractions, and the integral colors are well re- 
produced (Hernandez, Avila-Reese & Firmani 2001). 
The rotation curve decomposition is also reproduced, 
albeit with an excess of dark matter within the solar 
radius, Rq = 8.5 kpc. At Rq, the gas and stellar 
surface densities, as well as the B — K color at the 
present epoch, agree rather well with observational 
estimates (Fig. lb). 

The present-day global SFR of the MW model is 

2.9 M©yr _1 ; it increases slightly towards the past, 
attaining a maximum at 9.8 Gyr ago (z ss 1.3). The 
local surface SFH at Rq is shown in Figure la (thick 
solid line, in unities of M Q yr~ 1 pc~ 2 ). It starts to 
increase 1 1 Gyr ago (roughly when the disk forms at 
this radius), attains a maximum ~ 6 — 7 Gyr ago, 
and then decreases by a factor of two to the present 
epoch. By using a sample of field stars from the 
Hipparcos catalogue and comparing it with synthetic 
colour-magnitude diagrams, the SFR per unit of vol- 
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ume at Rq was infered by Bertelli & Nasi (2001). 
Hernandez et al. (2001) converted this SFR to a 
surface SFR (open circles in Fig. la). The agree- 
ment between observations and model is encourag- 
ing. The introduction of a soft core in the CDM halo 
(to fit better the observed rotation curve) does not 
influence significantly on the predicted SFR history 
(dashed line in Fig. la). We notice that the obser- 
vational inference describes only the general trend 
of the SFH over the total life of the system, and 
not its detailed shape, which may be discontinuous. 
The model results were smoothed in time to facili- 
tate the comparison. Interactions with satellites can 
introduce bursts of SF (e.g., Rocha-Pinto, Maciel & 
Scalo 2000; Kauffmann, Chariot & Balogh 2000). 

3.2. Integral disk population SF history 

Figure 2 reproduces a compilation of data for the 
integral SFR in the universe at different redshifts 
presented by Kravtsov & Yepes (2000). In this Fig., 
the integral SF history due to only (model) disks is 
also shown (thick dashed line). At low redshifts the 
observed SFR decreases more rapidly than models. 
Models give actually an upper limit for the disk con- 
tribution to the cosmic SF history, because they do 
not take into account neither the truncation of the 
gas infall in the massive galaxies due to large gas 
cooling times and/or external factors, nor the fact 
that at late epochs some of the disk accreting ma- 
terial might already be in form of stars (satellite or 
companion galaxies). At zZ 2, most of the data col- 
lected by Kravtsov & Yepes, as well as more recent 
estimates (e.g., Nandra et al. 2002; shaded rectan- 
gle), show that the SFR remains roughly constant, 
while the (model) disk contribution decreases. A re- 
cent re-analysis of the HST deep field by Lanzetta 
et al. (2002) shows that the SFR even increases sig- 
nificantly at redshifts larger than 2. The model disk 
galaxy population (ACDM cosmology) seems to be 
able to explain the observed cosmic SF history in- 
ferences until z ~ 2, modulo the strong fall since 
z w 0.5. At epochs earlier than z«2, the contribu- 
tion of disk galaxies to the cosmic SFR slightly and 
continuosly decreases. 

4. CONCLUSIONS 

The SF process is a key ingredient for models 
of galaxy formation and evolution. For disk galax- 
ies, stationary (non-bursting) and self-regulated SF 
seems to be a reasonable description. In this descrip- 
tion, the large-scale SF is triggered by gravitational 
instabilities of the gaseous disk (Toomre criterion) 
and self-regulated by an energy balance in the ISM 
perpendicular to the disk plane. The key factor in 
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Fig. 2. Dust corrected estimates of the cosmic SFR at 
different redshifts (points) from UV and Ha observations 
as collected by Kravtsov & Yepes (2000), and from X- 
ray observations of z ~ 3 LBGs. The data points have 
been converted to the ACDM cosmology. Dashed region 
is from the Kravtsov & Yepes gasdynamical simulations, 
while the dashed solid line is our model estimate of the 
contribution to the cosmic SF history by stars formed in 
the disks (the ACDM model was used). 

this energy balance is the ability of the turbulent 
ISM to dissipate the kinetic energy; the turbulent 
dissipation time determines the effectiveness of large- 
scale SF feedback. 

Our MHD, nonisothcrmal simulations of the ISM 
including realistic, stellar-like forcing, show that tur- 
bulent energy is dissipated rapidly, mostly at or near 
the sources, where shocks are common. The dissipa- 
tion time td is found to increase with the forcing 
scale If and depends only very weakly on e^ and S*. 
For properties of the ISM close to those at the solar 
neighborhood, td~ 15 — 20 Myr. Forced and decay- 
ing turbulent regimes coexist within the same flow. 
The rms velocity dispersion of the "residual" tur- 
bulence decays roughly inversely with distance from 
the injection sources. Turbulent motions produced 
near the disk plane will propagate only to distances 
of the order of the observed HI disk semi-thickness. 
This is consistent with the disk self-regulated SF 
scheme mentioned above, where the HI thickness re- 
sults from the balance between kinetic energy injec- 
tion and its dissipation rate. 

The application of this scheme to self-consistent 
models of disk galaxy formation and evolution within 
growing ACDM halos leads to rather succesful pre- 
dictions for the SFR, surface brightness, and color 
radial profiles of the disk galaxy population. The 
main drivers of the SF history are both the gas infall 
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rate determined by the cosmological mass aggrega- 
tion rate, and the disk surface density determined 
by the halo angular momentum. The SFR depends 
very weakly on mass. The SF history predicted for 
the solar neighborhood is in excellent agreement with 
observational inferences. The predicted integral disk 
SF history describes well the observed cosmic SF his- 
tory until z = 2 — 3. At earlier epochs, the disk pop- 
ulation SFR decreases while the observations show 
that the cosmic SFR remains constant or even in- 
creases, suggesting that at early epochs most of the 
UV and X-ray radiation is not produced by SF in 
disks. 

This work was supported by CONACyT grant 
33776-E to V.A. 
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